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PREFACE 


Final report of Asymmetrical Booster Ascent Guidance and Control System Design 
Studies performed under Contract NAS9-I3568 are contained in five separate 
volumes identified as follows: 

Volume I - Summary 

Volume II - SSFS Math Models - Ascent 

Volume III - Space Shuttle Vehicle SRB Actuator Failure Study 
Volume IV - Sampled Data Stability Analysis Program (SADSAP) 

Users Guide 

Volume V - Space Shuttle Powered Explicit Guidance 
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ABSTRACT 


This manual presents Boeing developed boost to orbit math models for the 
NASA/JSC Space Shuttle Functional Simulator. 

KEY WORDS 

Space Shuttle Vehicle 
Math Models 

SSFS (Space Shuttle Functional Simulator) 

Boost Dynamics 
Simulation Models 
Flight Dynamics 
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1.0 INTRODUCTION 

This manual presents the engineering equations and math models developed by 
the Boeing Aerospace Company for use in the Space Shuttle Functional Simulator 
(SSFS). These models were originally developed for NASA/JSC under Contract 
NAS9-12183, and continued under Contract NAS9-13568. This manual contains 
extensive revisions and additions to earlier documentation and it supersedes 
the previous math models document, Boeing Memorandum 5-2581 -HOU-1 02 dated 
4 October 1972. 

Section 2 contains documentation of all Boeing developed math models including 
several proposed models not yet incorporated into the SSFS. Included in 
section 3 are definitions of coordinate systems used by the SSFS models and 
coordinate transformations- 

Documentation of the flexible body math models is provided in section 4. 

These models have been incorporated in the SSFS and are in the checkout 
stage. 
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2.0 SSFS MATH MODELS 

This section contains environment math models for the SSFS computer program. 
Several model changes are contained within this documentation which have not 
yet been incorporated into the program. Subroutine ATTUDE has been deleted 
as a scheduled model and has been made a part of the flight control system 
subroutine. The control system now provides the call for attitude commands. 
Subroutine FLTSEQ (Flight Sequencing Program) has been included to simulate 
flight control system logic necessary to initiate staging signals. 

Each model is discussed under the following format: 

X.l Program Description 
X.2 Math Model 

X.3 Nomenclature (that used in math model) 

X.4 Input/Output Requirements' 
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2.1 ACCEL (Acceleration) 

2.1.1 Program Description 

This routine calculates acceleration in body coordinates for the accelerometer 
position APOS (X, Y, Z) on the vehicle. The component of acceleration returned 

is determined by the value of IFLAG. If IFLAG = 1 , the X component of accelera- 

) 

tion is returned. If IFLAG =2, the Y component of acceleration is returned. 

If IFLAG = 3, the Z component of acceleration is returned. All rotational 
and translational effects are included in the acceleration calculations. 

This routine is called by the control routine and is used to generate acceler- 
ometer signals for use by the flight control system. 
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2.1 .2 Math Model 


ENTER 


go to I FLAG 


calculate moment arms 


ACCEL = y j - + 

b> b> t 

X y y 

« b) - b) 

+ 0) 0) ii 

y X z X 

X Z 2 

:x 

+ 


eg 



Calculate X 
acceleration 


RETURN 


ACCEL • 

_ b) - 10 + o> (0 t 

z y X y X y X 


^cg 


Calculate Y 
acceleration 


RETURN 


firm = u)£ - wJi +(o(oJi 

fttULL y y X Z X 

_ b) - b) + u 

X z y z y z y 


+ 

eg 


Calculate Z 
acceleration 


RETURN 
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2.1.3 Nomenclature 


APOS^, APOS^, APOS 

X j ^ 

C6x> 



I FLAG 


Accelerometer location (X, Y, Z) 

Vehicle center of gravity (X, Y, Z) 

Moment arms (X, Y, Z) 

Vehicle center of gravity accel (X, Y, Z) 
Vehicle body angular rates (X; Y, Z) 

Vehicle body angular accelerations (X, Y, Z) 

Flag to specify component of acceleration 
to be returned 

I FLAG = 1 returns X component 
IFLA^ = 2 returns Y component 
IFLAG = 3 returns Z component 
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2.1.4 Input/Output 

This routine requires, as formal parameter input, accelerometer location and 
a code identifying the desired component of acceleration to be returned. Body 
angular rates, angular accelerations, translational accelerations, and C.G. 
locations must be input via common. Acceleration of the position denoted is 
output. 
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2.2 ACTVEH (6 DOF Equations of Motion) 

ACTVEH defines the motions of the center of gravity of the vehicle. For con- 
venience it is separated into five parts; 1) translation equations, 2) rotation 
equations, 3) euler angles, 4) initial position calculations, and 5) momentum 
transfer at staging. 

These first three equations should be solved at least once each second during 
powered flight. In the vicinity of environmental discontinuities more fre- 
quent solution is required; for instance, the vehicle can fly completely 
through a wind gust at maximum dynamic pressure within 0.1 second. Other 
discontinuities include*, staging, start of closed loop guidance, and engine 
or actuator failures. As a rule of thumb, the integration rate during 
transients can be 1/2 x rotational acceleration {in degrees/sec ). 

2.2.1 Translation Equations 

2. 2. 1.1 Program Description 

This model defines the linear accelerations of the rigid body. 

2.2.1 .2 Math Model 


EFy 


CD 

X 

u. 

1 

rFy 

P 

= [6] 

ZFw 

B 

Zfy 


SF^ 

^bJ 
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X 

1 


r 

9y 

T 

= [a] 

9y, 

h 

K 

1 

- « 


Xp = 9 j^ + zfy, /m 


/ ll ’‘“t ^ *( 


P ''P 




^P " ^Yp ^^Yp^"’ '^P ^ )t. ^Pjj 


xp. = ;,2 ildt . xp^ 

^p “ ftj 


^p ■ \ ^ ^P ■ /tj ^Pq ^P 


ill 


+ Zr 


2. 2. 1.3 Nomenclature 


zFy , zFw , zF, = sum of forces in the X, Y, Z body axis directions. 
^B ^B 


zFy , zFy , zF, = sum of forces in the X, Y, Z inertial plumbline axis 

'p 

directions. 

[B] = transformation matrix from body to inertial plumbline. 

zF = aero forces + thrust forces + RCS forces + engine deflection 
forces + slosh forces. 

9y » 9v » 9? “ gravitational acceleration components in inertial polar- 

,^I equatorial axis directions. 

9y » 9 y > 9? = gravitational acceleration components in inertial 

T ^P ^P plumbline axis directions. 

[a] = transformation matrix from inertial polar - equatorial to 

inertial plumbline. 

Xp. Yp, Zp = accelerations in inertial plumbline axis directions 
m = total vehicle mass 
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2. 2. 1.4 Input/Output 

The translation equations require as inputs: 

Aerodynamic forces 
Thrust forces 
RCS forces 

Engine deflection forces 

Slosh forces 

[a] and [g] matrices 

Gravitational acceleration components 

Vehicle mass 

* « • 

Initial conditions on Xp» ^p> ^p> 

The outputs from the translation equations are: 

zFv » zFw , ZF^ and the inertial plumbline position, velocity 
and acceleration components. 

The translation equations require the presence of subroutines: RCS, 
THRUST, AER0,TVC, SLOSH AND GRAVITY, 
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2.2.2 Rotational Equations 

2. 2. 2.1 Program Description 

This model defines the angular accelerations of the rigid body assuming 
that the center of mass lies approximately in the X-Z plane (ly^ = “ 0)- 

2. 2. 2. 2 Math Model 

^ **ZZ ■ 'xx* * ’xz 

p = (a ^ b I^)/C 

r (a + b ^yty)/^ 

a = zMy ^XZ 

b = zM, + pq - lyy) - qr 

D 

c = I I - I ^ 

^XX ^ZZ *xz 


2. 2, 2. 3 Nomenclature 


q = Angular acceleration about the Y body axis 
p = Angular acceleration about the X body axis 
r = Angular acceleration about the Z body axis 

^XX’ ^YY’ ^ZZ " inertia about X, Y, Z body axis respectively. 
Ixz = X “ Z Cross product moments of inertia 
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p, q, r = Integral of p, q, r (Body rates) 


iMy , IM, 


2 . 2 . 2. 4 
Inputs : 


Outputs: 


, 1M7 = Sum of moments about X, Y, Z body axes 

B 

= Aero moments + thrust moments + RCS moments 
+ engine deflection moments + slosh moments 


Input/Output 

From THRUST, AERO, RCS, TVC and SLOSH 

Moments (about body axes) due to aerodynamics, main propulsion, 
reaction control, engine accelerations and slosh. 

P, Q, R, P Q, R to IMU, Aero and Euler Angles 
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2.2.3 Euler Angles 

2. 2. 3.1 Program Description 

This model defines the rate of change of the euler angles describing the 

attitude of the vehicle in inertial space. 

2. 2.3.2 Math Model 

0 = (q cos (j) - r sin (j))/cos 
^ = q sin <|i + r cos 
^ = p = tan (p (q cos <)) + r sin(j)) 

2. 2. 3. 3 Nomenclature 

e, ^ = Euler angle rates (1st, 2nd, and 3rd rotations, respectively) 
e, = Integral ofe, i, ^ 

2. 2. 3. 4 Input/Output 

Inputs p, q, r from rotation equations 
Outputs e, to[B] and to IMU 
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2.2.4 


Initial Position Math Model 


2. 2. 4.1 Program Description 

This program calculates the difference between geodetic and geocentric 
latitude and uses it to calculate the initial state vector. This calcula- 
tion needs to be done once each time either the launch azimuth or latitude 
is changed. 

2. 2. 4. 2 Math Model 

i|^l^ = Arctan ^(1-f)^ tan j 
Re 

R. = /■ o N ^0 

Vl-f (2-f) cos'^ 


B 

Rv 


= 'J'L " '('l 

= cos e 
= Rj^ sin 6 sin 

= -R|^ sin 6 cos A^ 

= 0 

= a> cos A^ cos 
= u R|^ sin A^ cos 
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2. 2. 4. 3 Nomenclature 



= rotation rate of earth 
= launch geodetic latitude 
= earth flattening constant 

= earth equatorial radius 

= launch azimuth 
= launch geocentric latitude 
= magnitude of initial position vector 
= difference between geodetic & geocentric latitude 

= initial position vector in platform coordinate 

•= initial velocity vector in platform coordinate 

= altitude of vehicle CG above Fischer ellipse 


2. 2. 4. 4 Input/Output 

The constants needed by this model are w, f» R^, and A^. The 
output of the program is R^^, Ry> R 2 > Vy* ^Z* 


14 



2.2.5 


Staging Momentum Transfer 


2. 2. 5.1 Program Description 

This model accomplishes momentum transfer to the orbiter at the time 
of staging. 


2. 2. 5. 2 Math Model 


^0 


he 

r ^ 

ACGy 



= 

he 

. [BJ 

ACGy 


_Zo_ 


he 

^ACG^ . 


'V 


\c ' 

. [b] 

‘r^C-'CCy - 



he 

- '"Lc‘“x 

1-^0- 


L^lc- 

.j 



2. 2. 5. 3 Nomenclature 


Xq» Yq, Zq 

\c’ 


^ 0 ’ ^ 0 ’ ^0 


he \c’ ^LC 


[b] 


LC 


Orbiter inertial positions immediately after separation 

Launch configuration inertial position immediately 
before separation 

Orbiter inertial linear velocity components immediately 
after separation 

Launch configuration linear velocity immediately before 
separation 

Body to inertial transformation matrix 

Launch configuration angular rotation about X body axis 
immediately before separation 


Launch configuration angular rotation about Y body 
axis immediately before separation 

Launch configuration angular rotation about Z body axis 
immediately before separation 
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Xc’ ''“lc’ 

Center-of-gravity location in body coordinates for 
launch configuration immediately before separation 

’‘“o’ ’'“o’ ^“o 

Center-of-gravity location in body coordinates 
orbiter immediately after separation 

for 

‘“x 

pre-separation 


*CGy 1 

IVr “ Xr 1 pre-separation 



^CG^C ' 



2, 2. 5. 4 Input/Output 

Inputs: Vehicle Pre-separation C.G. Locations 

Vehicle Pre-separation Inertial Position 
Vehicle Pre-separation Velocity Components 
Vehicle Pre-separation Angular Rates 

from MASPRO and ACTVEH (translational and rotational sections) 
Outputs; Orbiter inertial position and velocity components 
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2,3 AERO (Aerodynamics) 

2.3.1 Program Description 

This model calculates and sums aerodynamic forces and moments for the 
vehicle. In addition this model calculates the latitude and longitude 
of the vehicle, flight path angle, mach number, dynamic pressure, angle- 
of-attack and sideslip angle, and the contribution to velocity due to 
wind speed and direction. 


2.3.2 Math Model 
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Vy = table lookup ~f (altitude) 
= table lookup ~f(al titude) 



1 
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Obtain from ATM0S 

1) p = f (altitude) 

2) a = f (altitude) 

3) P = f (altitude) 

H - Vj/a 
q = 1/2 P Vg2 

The following aerodynamic coefficients are looked up in tables for lift- 


off. 

boost to SRB 

separation. 

or orbiter 

+ ET 

to 

MECO. 

< 



^0 

= f 

(M) 

CL 

f 

(M) 

S' 

= f 

(M) 

r 

f 

(M) 

a 

= f 

(M, a) 

S' 

f 

(M) 

S' 

= f 

(M) 


f 

(M) 

\ 

= f 

(M) 

S' 

f 

(M) 


= f 

(M) 


f 

(M) 


= f 

(M, a) 


f 

(M) 


= f 

(M) 




a 



e 









Cy 

= f 

(M) 











"0 

Criteria for the selection of which group of tables to use should be 
as follows; 

1) If M < 0.6, use liftoff aero data tables. 

2) ■ If M > 0.6 and prior to SRB separation, use orbiter + ET + SRB 

aero data tables. 

3) If time is past SRB separation time, use orbiter + ET aero 
data tables. 
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Compute aerodynamic forces and moments 




F. = qS (C„ + C„ a) 

Ax Xq 

q S b 


F. = qS Cy 6 + 2 Vp Cy P + q S Cy 5r 
Ay RXg T 6r 

fA, = <’5 a ) + q S C^_ Se 


"0 a 


‘6e 


"''x = \ <^CG - - fA, (’CG - \> <'l, ^ " 

q S b^ 

C, sate, ar)a 2V„ (C, P + C, R) 


■5a 


6r 


\ "I 5 ^ “> * '’‘cG - ’'a„> ■ ''a^ <^CG - \> 

J 

4 q S 6e a Qa/2V,^B) 

”a. ' ^ ® ■ \ <’‘CG ■ ’'a^> ■" \ <''CG ■ 

q S - 

+ TV^ C P-qSb (C 6a + C 6r) 

Rx "p "6a "6r 


4 ' 

2.3.3 Nomenclature 

Xp. Yp, Zp = vehicle position in inertial polar-equatorial coordinates 

Xp, Yp, Zp * vehicle position in inertial plumbline coordinates 

[a] = transformation matrix from inertial polar-equatorial to 

plumbline coordinates 

Ay = latitude of present position of vehicle 

= East longitude of present position of vehicle corrected 
for earth's rotation 

t^ = time of launch (from epoch) i 

= radians to degrees conversion constant 
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w 


e 


t 

R 


V 


'^EARTH^ 

'^EARTHy 



= angular rate of earth 
= elapsed time from liftoff 

= distance from the center of the earth to the vehicle 

= components of earth's velocity in inertial polar- 
equatorial coordinates 

= components of vehicle relative velocity in plumbline coordinates 


Xp» Yp, Zp = components of vehicle velocity in plumbline coordinates 



= total vehicle relative velocity in plumbline coordinates 
= components of relative velocity in local vertical coordinates 


fol = transformation matrix from local vertical to inertial polar- 
equatorial coordinates 



= vehicle flight path angle with respect to local horizontal 
= horizontal wind speed in local vertical coordinates 
= wind azimuth (North * 0°) 

= components of wind velocity in plumbline coordinates 


= vehicle velocity with respect to air in body coordinates 



= transformation matrix from body to plumbline coordinates 


Vg = total vehicle velocity with respect to air in body coordinates 
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vehicle angle of attack 
vehicle sideslip angle 
local air mass density 
local speed of sound 
local air pressure 
Mach number 
dynamic pressure 

components of aerodynamic force in body coordinates 

vehicle aerodynamic reference area 
vehicle mean aerodynamic chord 
vehicle reference span 
aileron deflection 
elevator deflection 

rudder deflection 
vehicle roll rate 
vehicle pitch rate 
vehicle yaw rate 

Aerodynamic moments about the X, Y, and Z body axes 
respectively 



2.3.4 Input/Output 
Input from routines: 


Xp, Vp, Zp 

vehicle position in inertial plumbline 
coordinates from EOM 

Xp. Yp, Zp 

vehicle velocity in inertial plumbline 
coordinates from EOM 


wind velocity from tables 

AZ„ 

wind azimuth from tables 

p > a» p 

current air density, speed of sound and air 
pressure from ATMOS 

^CG’ ^CG’ ^CG 

current location of vehicle center of gravity 
from MASPRO 

t 

elapsed time from liftoff from flight sequencer 

Rv 

distance from center of the earth to the vehicle 
from E0M 


aerodynamic control surface deflections from 
flight software commands 

P, Q, R 

vehicle roll, pitch and yaw rates from E0M 

All aerodynamic 

coefficients are input from tables. 

Input from cards 

for initialization: 


time of launch (from epoch) 

Cl 

radians to degrees conversion constant 


angular rate of earth 

s 

vehicle aerodynamic reference area 

C 

vehicle mean aerodynamic chord 

b 

vehicle reference span 

Mr Mr Mr 

aerodynamic reference location in body 
coordinates 
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Output to routines: 




MX tty 


Output to printer: 
X.. 


V 


X 


LV 


Y 

AZ 


W 


V 

a 


B 


P 

a 


P 

M 

q 



current air pressure to THRUST 
components of aerodynamic forces to the E0M 

moments due to aerodynamic forces to E0M 


latitude of vehicle's position 
longitude of vehicle's position 
rate of climb 

flight path angle 
wind speed 
wind azimuth 

vehicle velocity with respect to air 

angle of attack 

angle of sideslip 

local air mass density 

local speed of sound 

local air pressure 

Mach number 

dynamic pressure 

aerodynamic control surface deflections 
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2.4 ATMOS (Atmosphere) 

2.4.1 Program Description 

This program calculates the speed of sound, pressure and air density 
from an altitude input, 

2.4.2 Math Model 

Use the Cape Kennedy Reference Atmosphere (TM-X-53872, PARAGRAPH 14.7 
MSFC "COMPUTER SUBROUTINE PRA-63") as specified for SSV design studies 


2.4.3 Input/Output 

The altitude above the mean earth surface must be supplied to the 
model which returns the speed of sound, pressure, and atmospheric 
density. 
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2.5 BLC (Baseline Control System) 

2.5.1 Program Description 

This model issues commands to the engine gimbals (to subroutine TVC via 
CMDFIL) such that the actual vehicle attitude is made to follow the 
attitude prescribed by the guidance model. This documentation represents 
the implementation of the RI system as described in the July 73 Space 
Shuttle Guidance and Control Data Book, with the addition of roll rate 
crossfeed into the Y accelerometer channel. With proper input of gain 
tables, this routine may be used both for first and second stage control. 
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2.5.2 Math Model 


( ENTER ) 


IF (tj 


> t > t„^) 


then 


AOP 

(1-8) = 

AOY 

(1-8) = 

RP 

(1-7) - 

RR 

(1-2) = 

RY 

(1-7) = 


= J 


RAMP TO 
UPPER 
STAGE 
VALUES 


PRE-SEPARATION 

RAMP MIXING 
LOGIC COEFS. 

TO UPPER 

STAGE 

VALUES 


^ ^ 

0^- = Interpolated from tables UNI 
AND FROM GUIDE THEREAFTER 

'I'^T = Interpolated from tables UNI 
AND FROM GUIDE THEREAFTER 

d>^T = Interpolated from tables UNI 
AND FROM GUIDE THEREAFTER 

= Interpolated from tables 

•IL GUIDANCE START - 
■IL GUIDANCE START - 
‘IL GUIDANCE START - 


L- f 

-KAt * , . /T -KAtx . , 

(|)Ci = e * 'i-Qj + (1-e ) + ({.j 

sj (test pass) 


DETERMINE 

COMMAND 

SIGNALS 


FILTER ROLL COMMAND 
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INTERPOLATE CONTROL 
GAINS AND ENGINE 
TRIM COMMANDS AS A 
FUNCTION OF TIME 
OR VEHICLE RELATIVE 
VELOCITY 


CALCULATE INERTIAL 
ATTITUDE ERRORS 


h 


TRANSFORM TO 
BODY ERRORS 


CALL FLIGHT SEQUENCING LOGIC 
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SCALE AND FILTER 
ROLL RATE SIGNAL 



SUM AND SCALE 
YAW ATTITUDE, 
RATE ACCEL AND 
ROLL RATE CROSS 
FEED SIGNALS 


SUM AND SCALE 
ATTITUDE AND 
RATE SIGNALS 
FOR SRB 


SUM AND SCALE 
ATTITUDE AND 
RATE SIGNALS 
FOR ORBITER 
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SCALE AND LIMIT 
ROLL ATTITUDE 
SIGNALS 


SUM ATTITUDE 
AND RATE SIGNALS 


SCALE SIGNAL 
FOR SRB 


SCALE SIGNAL 
FOR ORB I TER 


PRODUCE LIMITED 
ORBITER SIGNAL 
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"^OPl 

^oP2 

^oP3 

■^oYI 

^oY2 

^oY3 


= AOP(l) * PORB + RP(1) * R0RB5 + TRIM^ 

** A0P{2) * PORB + RP(2) * RORB5 + TRIM^ 

= A0P(3) * PORB + RP(3) * RORB5 + TRIH^ 

= AOY(T) * YORB + RY(1) * RORB 

= A0Y(2) * YORB + RY(2) * RORB 

= A0Y(3) * YORB + RY(3) * RORB 



MIX SIGNALS TO GENERATE 
INDIVIDUAL ORBITER 
ENGINE COMMANDS 
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G 


= A0P{4) * PSRB + RP(4) * RSRM + TRIM^ 

= A0P{5) * PSRB + RP{5) * RSRM + TRIM^ 

= A0Y(4) ( YSRM + RY(4) * RSRM 

= A0Y(5) ( YSRM + RY(5) * RSRM 


MIX SIGNALS TO GENERATE 
INDIVIDUAL SRB ENGINE 
COMMANDS ASSUMING ACTUATOR 
ORIENTATION IS 0“ - 90° 


DELP4 * COS 45° - DELY4 * SIN 45‘ 


DELP5 * COS 45° - DELY5 * SIN 45' 


DELP4 * SIN 45° + DELY4 * COS 45' 


= DELP5 * SIN 45° + DELY5 * COS 45' 


TRANSFORM SIGNALS TO 
GENERATE SRB ENGINE 
COMMANDS FOR 45° - 135" 
ACTUATOR ORIENTATION 
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GENERATE RUDDER 
COMMAND 
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2.5.3 

AOP(i) 

AOY(i) 

K 

KAQ 

KAR 

KT 

KY 

KZ 

K(|) 

•■o 

PORB 

PSRB 

Qk 

Qo 

Q4 

Q 5 


RA 

RAp 

RORB 

R0RB5 


Nomenclature 

Mixing Logic Coefficients for Pitch Actuators Used for 
Pitch Control, i = 1-8 

Mixing Logic Coefficients for Yaw Actuators Used for Yaw 

Control , i = 1-8 

Roll Command Filter Constant 

Pitch Attitude Error Gain 

Yaw Attitude Gain 

Rate Gain 

Lateral Accelerometer Gain 
Normal Accelerometer Gain 
Roll Attitude Gain 
Orbiter Roll Rate Gyro Signal 
Pitch Signal to Orbiter Engines 
Pitch Signal to SRB Engines 
Summed SRB Pitch Rate Signal 
Orbiter Pitch Rate Gyro Signal 
Left SRB Pitch Rate Gyro Signal 
Right SRB Pitch Rate Gyro Signal 
Orbiter Yaw Rate Gyro Signal 
Left SRB Yaw Rate Gyro Signal 
Right SRB Yaw Rate Gyro Signal 
• Roll Rate Crossfeed Signal 
Filtered Roll Rate Crossfeed Signal 
Roll Signal to Orbiter Engines 

Roll Signal to Orbiter Engines Limited to ±5 degrees 
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RP(i) 


RR(i) 

RSRB 

RUDDER 

RY(i) 


TRIMq 

TRIM3 

Y 

h 

^FO 

YORB 

YSRBL 

YSRBR 

Z 

Z„ 


^FO 

^CI 

‘^CI 

***CIF 


Mixing Logic Coefficients for Pitch Actuators Used for 
Roll Control, i = 1-7 

Mixing Logic Coefficients for Aero Actuators Used for 
Roll Control, i = 1-2 
Roll Signal to SRB Engines 
Deflection Command to Rudder 

Mixing Logic Coefficients for Yaw Actuators Used for 

Roll Control, i = 1-7 

Trim Signal to Orboter Pitch Signal 

Trim Signal to SRB Pitch Signal 

Summed Lateral Accelerometer Signal 

Aft Lateral Accelerometer Signal 

Filtered and Scaled Lateral Accelerometer Signal 

Forward Lateral Accelerometer Signal 

Yaw Signal to Orbiter Engines 

Yaw Signal to Left SRB Engine 

Yaw Signal to Right SRB Engine 

Summed Normal Accelerometer 

Aft Normal Accelerometer Signal 

Normal Accelerometer Command 

Filtered and Scaled Normal Accelerometer Signal 

Forward Normal Accelerometer Signal 

Inertial Flight Path Angle Rate Command 

Inertial Roll Attitude Command 

Filtered Inertial Roll Attitude Command 
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0 

e 

0 

0 

0 

0 

0 

f 

V 
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V 
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6 

6 

6 

6 

6 
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eB 

el 

I 

K 

S 

K 

Cl 

eB 

el 

I 

K 

S 

CB 

eB 

el 

I 

CB 

K 

S 

oPl 

oP2 

oP3 

SR5 

SR6 

oYl 

oY2 


Body Roll Error 

Inertial Roll Error 

Inertial Roll Attitude 

Scaled Roll Attitude Signal 

Summed Roll Attitude and Rate Signal 

Roll Attitude Signal Rate of Change 

Inertial Pitch Attitude Command 

Body Pitch Attitude Error 

Inertial Attitude Error 

Inertial Pitch Attitude 

Scaled Pitch Attitude Signal 

Summed Pitch Attitude and Rate Signal 

Body Pitch Rate Command 

Body Yaw Attitude Error 

Inertial Yaw Attitude Error 

Inertial Yaw Attitude 

Body Yaw Rate Command 

Scaled Yaw Attitude Signal 

Summed Yaw Attitude and Rate Signal 

Deflection Command to Pitch Actuator of Orbiter Engine #1 

Deflection Command to Pitch Actuator of Orbiter Engine #2 

Deflection Command to Pitch Actuator of Orbiter Engine #3 

Deflection Command to Right Actuator of Left SRB 

Deflection Command to Right Actuator of Right SRB 

Deflection Command to Yaw Actuator of Orbiter Engine #1 

Deflection Command to Yaw Actuator of Orbiter Engine #2 
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^SL5 

*^516 

^OMSp^ 

*0MSp2 

^OMSy-, 

^OMSyg 

At 


Deflection Command to Yaw Actuator of Orbiter Engine #3 
Deflection Command to Left Actuator of Left SRB 
Deflection Command to Left Actuator of Right SRB 
Deflection Command to Pitch Actuator of OMS Engine #1 

Deflection Command to Pitch Actuator of OMS Engine #2 

Deflection Command to Yaw Actuator of OMS Engine #1 

Deflection Command to Yaw Actuator of OMS Engine #2 

Computation Frequency for Control Subroutine 
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2.5.4 Input/Output 

Inputs to this model are inertial attitude angles, body rotational rates, 

Y and Z translational acceleration, prestored engine trim deflection 
commands, prestored acceleration commands, and attitude commands from guidance. 
Outputs from this model are engine gimbal and aero surface deflection commands 
output via the CMDFIL routine. 


40 



2.6 CGAINS {Control Gains Equations) 

2.6.1 Program Description 

The CGAINS program is used to calculate the control gains necessary for a 
desired type of control during the Shuttle boost. There are several options 
for the control gains that are calculated: load minimum, drift minimum, or 

attitude control for the pitch and yaw gains; and thrust vector control or 
aero control for the roll gains. The following model presents the equations 
necessary to calculate these gains. 


2.6,2 Math Model 

The following quantities must be calculated each time the control gains 
are needed. The symbols used in these equations are defined in Table I. 
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Pitch and Yaw Control Gains for Load Minimum Option 
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Pitch and Yaw Control Gains for Drift Minimum Option 
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Thrust Vector Control 


Cpr = - [ S (Fx * RPi * + F * RY. 

2 


u. 


Or - 
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Aero Control 


2 c tu 
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C2r = -qSB [C^_ RRi -1 RR2] / Ix 




u 

a« = X 
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2.6,3 


Nomenclature 


Variable 

AOP. 

AOY^. 


^ly 


Definition 

Mixing Logic Coefficients for Pitch Actuators Used 
for Pitch Control, i = 1-8 

Mixing Logic Coefficients for Yaw Actuators Used 

for Yaw Control, i = 1-8 

Pitch Attitude Gain 

Yaw Attitude Gain 

Pitch Attitude Rate Gain 

Yaw Attitude Rate Gain 


b Wing Span 

C Coefficients for Calculation of Gain 


c 

C. 


Cm 


Cm, 


Cn. 

V 


Mean Aerodynamic Chord 

Aero Rolling Moment Coefficient Due to Aileron 
Deflection 

Aero Rolling Moment Coefficient Due to Rudder 
Deflection 

Aero Pitching Moment Coefficient Due to Angle-of -Attack 

Aero Pitching Moment Coefficient Due to Elevator 
Deflection 

Aero Yawing Moment Coefficient Due to Sideslip-Angle 
Aero Yawing Moment Coefficient Due to Rudder Deflection 

Aero Side Force Coefficient Due to Sideslip-Angle 

Aero Side Force Coefficient Due to Rudder Deflection 
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Variable 



®2p 

92y 

I 

K 

1 

m 

N' 

q 




s 

^ACCEL 

^CG’ ^CG’ ^CG 

Xr, Yr, 2r 



Def i ni ti on 

Aero Normal Force Coefficient Due to Angle-of -Attack 
Aero Normal Force Coefficient Due to Elevator Deflection 

Thrust in X direction from engine i 

Pitch Acceleration Gain 
Yaw Acceleration Gain 
Moments of Inertia 

Coefficients for Calculation of Gains 

Moment Arms 

Mass 

Partial of Normal Force 
Dynamic Pressure 

Mixing Logic Coefficients for Pitch Actuators Used 
for Roll Control, i = 1-7 

Mixing Logic Coefficients for Aero Actuators Used 
for Roll Control, i = 1-2 

Mixing Logic Coefficients for Yaw Actuators Used 
for Rol 1 Control , i = 1 -7 
Aero Reference Area 
X Accelerometer Location 

Location of Center of Gravity 

Location of Aero Reference Point 

Engine Locations, i = 1-6 
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Variable 





01 


Definition 
Temporary Variable 

Damping Ratio 
Natural Frequency 


2.6.4 Input/Output 

The control gains routine calculates the attitude gains, attitude rate gains 
and accelerometer gains for each of the desired conditions which have been 
previously . menti oned . These gains are also calculated for a range of natural 
frequencies (w). These results are output on a scratch tape or FASTRAND 
file for processing by a plotting program. 
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2.7 FLTSEQ (Flight Sequencing) 

2.7.1 Program Description 

Although not currently active, this program is called by the Flight Control 
System (BLC) and is used to initiate staging and provide flight sequencing 
logic as required. 

The switching points between stages were defined by the following events: 

1) SRB separation and guidance initiation; 2) main engines shutdown and RCS 
activation; 3) external tank jettison; 4) QMS burn initiation and RCS de- 
activation; and 5) OMS shutdown and RCS actuation at orbit insertion. 


2.7.2 Math Model 


BOOST TO ORBIT 
FLIGHT CONTROL LOGIC 
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ENGINE 

SHUTDOWN 




ACTIVATE RCS 




SWITCH TO "COAST 
STAGE BEFORE 
ET JETTISON" 




BOOST 


(FROM 

GUIDANCE) 


(TO ENVIRON) 


SET SEPARATION 
TIMER 





1 

' (TO ENVIRON) 

J 

)Si. 

SEt OMS 
START TIMER 




COMMAND ET 
SEPARATION 




SWITCH TO 
COAST STAGE AFTER 
ET JETTISON 










OMS BURN 



(FROM GUIDANCE) 


(TO ENVIRON) 



ACTIVATE RCS 


SWITCH TO 
"ORBIT COAST" 




2.8 GUIDE (Guidance Model) 

Due to recent revisions and updates to the guidance equations, these models 
are not shown here but are presented in Volume V of this report. 
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2.9 MASPRO (Mass Properties) 

2.9.1 Program Description 

This model provides the mass properties, which consist of center of gravity 
travel, moments of inertia as a function of weight and total mass calculation. 

2.9.2 Math Model 
W = W - 

At Staging W = W - W (I) 

CG = table lookup F(W) 

A 

CG^ = table lookup F(W) 

CG^ = table lookup F(W) 

I = table lookup F{W) 

A A 

I = table lookup F(W) 

^2Z " "tsble lookup F(W) 

I = table lookup F(W) 

2.9.3 Nomenclature 

W = total vehicle weight 
W(l) = solid rocket booster weight 
W(2) = external tank weight 
ft = weight flow rate input from THRUST 
CG = X center of gravity location from vehicle reference 

CGy = Y center of gravity location from vehicle reference 

■ CG^ = Z center of gravity location from vehicle reference 

I = Mass moment of inertia about X axis 

A A 

I = Mass moment of inertia about Y axis 

I = X - Z cross product of inertia 

zz 
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2.9.4 Input/Output 
Input: W from THRUST 

Output: W, CG's and I’s. 
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2.10 MAXMIN (Max and Min Parameter Dump) 

2.10.1 Program Description 

This program is used to collect the extreme values of certain flight parameters. 
At the stop time scheduled for MAXMIN the collected values are printed. This 
routine may be scheduled as often as desired, as long as the calls do not 
overlap. 

2.10.2 Math Model 
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2.10.3 Nomenclature 


V(i) 

MAX(i) 

MIN(i) 

ABS(i) 

TMAX(i) 

TMIN(i) 

TABS(i) 


Parameter value to be tested for magnitude 

Maximum value of parameter i collected 

Minimum value of parameter i collected 

Absolute maximum value of parameter i collected 

Time of occurrence of maximum 

Time of occurrence of minimum 

Time of occurrence of absolute maximum 



2.11 ORBITl (Orbit Parameters) 

2.11.1 Program Description 

This math model calculates the parameters of the trajectory achieved 
at insertion from knowledge of the state vector in polar equatorial 
coordinates at the time of orbiter engine shutdown. The program 
calculates node, inclination angle, orbit phase angle, eccentricity, 
orbit parameter, true anomaly, apogee altitude, and perigee altitude. 

2.11.2 Math Model 


n = Arctan 


~^YPE 

■^XPE 



+ VzPE ^YPeN 

^ ''ZPE ^XPE j 


I, = Arctan 


( 


*ZPE 

-Vxp£ Sinn + Vyp£ cos 






cosn ^ ^YPE 

V - Arctan^ cos ^(-Vxp£ sinn + Vypp cosn) + sin c 


'YPE 


7PE 




R = A / R^p£^ + RypE^ + R2pe^ 
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2 


jR X V| 


p = 

K 


K 

SMA = 


e = 

" SMA 

FP = 

R ’ V 

^rV - IR-V 


TA s Arctan 


Rp = SHA (1-E) - R^ 


R^ = 2 • SMA - Rp - 2 Rg 


2.11.3 Nomenclature 

'^XPE’ '^YPE’ '^ZPE ~ polar-inertial coordinates 

Rxpe» ^yPE’ ^ZPE "" Position in polar-inertial coordinates 
n = Longitude of the ascending node 
C = Inclination angle 
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e = Eccentricity 

K = Gravitational constant = 1.407654 x 10^^ ft^/sec^ 

V = Orbit phase angle (angle between equator and perigee) 

P = Semi-latus rectum 

SMA = Semi -major axis 

FP = Tangent of flight path angle 

TA = True anomaly 

Rp = Altitude at perigee 

= Altitude at apogee 

R^ = Mean radius of the earth 
e 

2.11.4 Input/Output 

The ORBITI math model requires the vehicle state vector in platform 
equatorial coordinates as input. As output the model prints out 
ascending node, inclination angle, orbit phase angle, eccentricity, 
orbit parameter, true anomaly, apogee altitude, and perigee altitude. 
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2.12 


ORBITR (3D EQUATIONS OF MOTION)' 

2.12.1 Program Description 

This Math Model calculates vehicle acceleration from active guidance 
commands and integrates to get velocity and position in platform 
coordinates. Polar-Equatorial and Local Vertical coordinates 
systems are erected to calculate latitude, longitude, and flight path 
angle. Logic is included for acceleration limiting, integration 
cycle time rectification, and velocity cutoff. 

2.12.2 Matli Model 



N 

E 

1=1 


Ty (I) 



_ P 

‘>AP 'TTS^ 

Apt' = Qpj Pj Q;^p + Ag^^^y 
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V 


V 


• it ^PT <lt + 

./t V. dt + R, 


'LV 


C«j'^ Rp' 


= Arcs! 


Rp(2) 


= Arctan ( ^ 


[6f [a 3^ Vp 


= Arcs in 


\,{ 1 ) 




LV' 


Integration Cycle Time Rectification 

Acceleration limiting will most likely occur between integration time 
points. Therefore the integration for this interval must be done in 
two parts . 

Velocity Cutoff 

9 

The program will be terminated either when |Vp| exceeds Z or when 
is less than whichever occurs first. 
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2 . 12.3 


Nomenclature 


\ 





^GRAV 










I 

M, 


= Thrust acceleration at full throttle 
= Number of engines 

= Maximum vacuum thrust of main engine 
= Current vehicle mass 
= Maximum allowed vehicle acceleration 
= Magnitude of acceleration due to thrust 
= Throttle setting 

= Total mass flow rate for all engines at full throttle 
= Acceleration command from guidance in platform coordinates 
= Unit vector acceleration command 

= Current total vehicle acceleration in platform coordinates 
= Acceleration due to gravity in platform coordinates 
= Velocity of vehicle in platform coordinates on last pass 
= Current vehicle velocity 

= Position of vehicle in platform coordinates on last pass 
= Current vehicle position 

= Position of vehicle in polar equatorial coordinates 

= Geodetic latitude of launch site 

= Longitude of launch site 
= Rotation rate of earth 
= Time of launch 
= Elasped time since launch 
= Mass of vehicle on last pass 
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A|^ = Launch azimuth 

Xyj = Present vehicle latitude 

= Temporary variable 

[а] = Transformation from polar equatorial to platform coordinates 

[б] = Transformation from local vertical to polar equatorial 

coordinates 

^ = Present vehicle longitude 

Y = Flight path angle 


^LV ~ local vertical coordinates 

2 = Target velocity in guidance coordinates 

MyQ = Mass of empty orbiter 


2.12.4 Input/Output 

The targeting program provides Z. The Active Guidance program 
provides ^p. The gravity program provides Agp^y. The resulting quan- 
titles calculated by this model are Xy» ^pj* ^LV* 

Flight software commands accepted by this model are Ap. The Active 
Guidance, Targeting, and Gravity programs must be present to provide 
inputs to the model. 
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2-13 ORBTAR (Boost Orbit Insertion Targeting Model) 

2.13.1 Program Description 

The Targeting program is used in the flight software to describe the 
orbit plane with respect to the launch pad. Position, velocity, and a 
unit vector normal to the orbit plane at perigee are, calculated from a 
knowledge of perigee and apogee altitudes, location of the launch pad, 
orbit inclination, and an orbit parameter. 

2."1T. 2 “Math Model ’ " " ' ‘ 



R = Rp 
Y = 0 


Z = unconstrained 
R = 0 
Y = 0 
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Sg = sin (x) cos(&)+ cos (A^) cos (x) sin (s) 

Cj^ = sin (A^) cos(x) 

Cg =Vl - Sg^' 



.. a = COS-(-Y)/Cg 

A = Arcsin (Cp) - Arcsin (a) 

Up (1) = -Sin (3) sin (a) 

Up (2) = cos (A) 

Up (3) = cos (3) sin (a) 


2.13.3 Nomenclature 

Rp = distance from earth center of mass to periapsis 
= distance from earth center of mass to apoapsis 
hp = altitude at periapsis 
h^ = altitude at apoapsis 
A = semi -major axis 
e = eccentricity 
Rg = radius of earth 
R = radial distance at insertion 
Y = cross-range distance at insertion 
Z = downrange distance at insertion 
R = radial rate at insertion 


69 



Y = lateral velocity at insertion 

Z = downrange velocity at insertion 
Sg, Cj^, Cg» Cp, a, A = temporary variables 

^ = latitude of launch pad 

M = universal gravitational constant 

3 — =-orbi-t-parameter 

= launch azimuth 

Y = orbit inclination angle 

Uq = unit vector normal to desired orbit plane in platform 
^ coordinates (see Section IV) 

2.13.4 Input/Output 

This model requires R^, hp, h^,A , B, A^, and Y as input and calculates 
* * * 

R, Y, Z, R, Y, Z, and Uq. The model needs to be called only once per 
simulation. 
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2.14 RCS (Reaction Control System) 

2.14.1 Program Description 

This program calculates the moments and linear accelerations applied 
to the vehicle due to RCS thrust commands from the flight software. 

2.14.2 Math Model 

Thrust commands are conditioned according to the following limitations: 

(1) A jet cannot be commanded to ignite unless the duration of 
ignition is some minimum value 

(2) No jet can be ignited continuously longer than Tg seconds 

(3) A jet cannot be commanded to ignite unless T^ seconds has 
elapsed since the previous ignition has ceased. 

The following shov-zs a logic diagram which can be used to Implement these 
limitations. 
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RCS LOGIC DIAGRAM 
FIGURE 2-1 
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2,14.3 Nomenclature 


i = engine index 

j = 1, 2, or 3 for x, y, or z respectively 
^MAX ” thrust achieved with continuous ignition 

^MIN “ average thrust for a computation cycle for minimum thrust duration 
TC^. (j) = commanded thrust 
T^(j) = realized thrust 

IM.(j) = indicator for continuous thrust for the previous pass 
' (set equal to zero for restart) 


IN.(j) = thrust enable flag > 

' (set equal to zero for restart) 

= check time for 

Thc " check time for Tg 

t = current time at entry to RCS program 




m;. = X T^' (1) 

P^(l) = ELX^. ^ XCG 

P.(2) = ELY^. - YCG (2) 

P.(3) = ELZ^ = ZCG 

P^- = engine position vector for the i jet; ELX^, ELY^. , ELZ^ _ 

X, Y, and Z locations, respectively, of the i^^ RCS jet cluster; 

XCG, YCG, ZCG = X, Y, and Z locations, respectively of the vehicle 
center of mass 
i ’= engine index 

— ^ th 

M^. = moment vector due to thrust from the i jet 
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T^. = thrust vector for the jet 
X = indicates vector cross product 

N 

MR = ^ 5 1 

MR = total effective moment from all RCS jets 
N = number of RCS jets 


N 



F = total effective linear force 

3 N 
''' "j?! 1 = 

A = reduction in mass of vehicle due to RCS fuel usage 
R = RCS jet efficiency constant 


( 3 ) 


(4) 


2.14.4 Input/Output 

The model requires ELX, ELY, ELZ, XCG, YCG, ZCG, T , t, and M as input. 
The flight software command accepted is T . The cross product, mass 
properties, and EOM subroutines must be present. The model provides MR, 
A, and F as output. 


74 



2,15 THRCMD (Throttle Command) 

2.15.1 Program Description 

This model senses vehicle longitudinal acceleration and issues throttle com- 
mands to the SSME’s if the "3g" acceleration limit is exceeded. The throttle 
command is calculated such that the limit is maintained. 

2.15.2 Math Model 








2.15.3 Nomenclature 

/\CCEL Vehicle acceleration along X body axis 

F Total thrust forces along X body axis 

X 

M Vehicle mass 

go Accel, of gravity 

F Allowable thrust to maintain 3g acceleration 

LIMIT 

Sum of SSME X axis thrusts 

SSME 

F Sum of SRB X Axis thrusts 

SRB 

THROT Throttle command to SSME's 

2.15.4 Input/Output 

This model requires as input X axis thrusts for the SRB and SSME 
addition the present throttle setting must be supplied. Outputs 
of identical throttle commands to the Space Shuttle main engines 


s. In 
consist 
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2.16 THRUST 

2.16.1 Program Description 

This model calculates forces and moments in body coordinates due to thrust 
forces from all engines given engine gimbal angles, throttle settings, 
and atmospheric pressure. 

2.16.2 Math Model 

Table lookup for SRM vacuum thrust, SRM mass flow rate, and power-on base 
drag. 

POBD = f (h) 

ti""esRM ■ 

T., = f (time<.R„) * (1 + UNBAL) 

''i (i = 4,5) 

^i (i = 4-5) " ^ ^^^'^^SRM^ * (1 + UNBAL) 


Calculate engine forces 


^i{i 

= P * 

= 1-7) 


- 


*»’a) 

®p = 

“2 * 

(tan 6p 
^i 

+ ai 

★ 

(tan 

6y )) 

0y = 

• 4 ( 

“2 

(tan 6w 

U 

■ “1 

★ 

(tan 

6p )) 

^i 


if i 

= 1-3, 6- 

■7 

“1 

= 0 

and tt 2 = 1 


i 

= 4,5 


“1 

= 1 

and og = .7071068 


(for i * 4,5; and Og are unique for diag, and SRB actuators) 

2 2 T ^/2 

TBX. = T. * (1 + q / + ey^) 

TBY^ = TBX^ * 0y 
TBZ. = TBX. * 0p 
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Sum Forces 


FBy = E TBX. + POBD 
^ i=l,7 

FB. = E TBY. 

^ 1=1.7 

FB, = E TBZ. 

^ 1=1,7 ’ 

Calculate moments 

MTXB. = TBZ. (ELY. - YC6) - TBY. (ELZ. - ZCG) 

MTYB. = TBX. (ELZ. - ZCG) - TBZ^ (ELX. - XCG)' 

MTZB. = TBY^ (ELX^ - XCG) - TBX^ (ELY^ - YCG) 

Sum moments 
MTBj^ = eMTXB^. 

MTBy = eMTYB. 

MTB^ = eMTZB. 

2.16.3 Nomenclature 

1 = engine index; 1,2,3 = orbiter engines, 

4,5 = solid engines 
6,7 = QMS engines 

POBD — Power-on base drag 

h Altitude 

UNBAL = Individual SRM thrust unbalance 
Ty = Vacuum thrust 

|5t^ = Mass flow rate 

T. = Altitude thrust 

1 
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Throttle setting 





Y 

TBX 


i 


TBY. 


TBZ. 



Engine exit area 
Atmospheric pressure 
Engine pitch gimbal angle 

Engine yaw gimbal angle 

Engine force in X body coordinate 
Engine force in Y body coordinate 
Engine force in Z body coordinate 
Total thrust force in X body coordinate 
Total thrust force in Z body coordinate 


MTXB^, MTYB-, MTZB^. = X, Y, and Z components* respectively, of moments 
due to engine number i, 

ELX^ , ELY^ , ELZ^ = X, Y, and Z components, respectively, of engine locations 
in body coordinates. 

XCG, YCG, ZCG = X, Y, and Z components of location of center of mass of 
the vehicle. 

MTBj^, MTBy, and MTB^ = X, Y, and Z components, respectively of total 
moments due to engine thrust 


2.16.4 Input/Output 

The parameters which must be supplied to the model as input are EA^ , ELX^ , 


ELY., ELZ 



, XCG, YCG, and ZCG. The resulting 

T - 
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quantities calculated by the model are FTB and MTB. The flight software 

command accepted by the model is ■ The outputs of this model are accepted 

i 

by the equations of motion. 
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2.17 TSHAPE (Trajectory Shaping) 


2.17.1 Program Description 

The trajectory and control parameters which must be calculated to accom- 
plish trajectory shaping are a^, desired angle-of-attack; desired 
pitch attitude angle; 6^, desired engine deflection angle; and 
desired body sensed acceleration. The values of e^, 6^. and are 
dependent on The boost flight is divided into three phases: vertical 

rise, tilt maneuver, and alpha policy. For each of these flight phases, 


0 ^ is calculated differently. 

2.17.2 

Math Model 

2.17.3 

General Calculations 


= C^a * 180 
z r L 

\ = q S 
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■ 

»z 

= F 

ax 

^CG 


*19 

II 

X 

o 

cr> 

-*o 


^19 

" ^CG 

, ■> 'i 
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L 


B 


^CGB 


“^CGO 







= -tan 


^CG - 1 

_*CG ■ *B , 


= tan 




Angle-of-Attack Calculations 


Vertical Rise 



®D 



~ Fq cos 6p = Fg cos 6g 
_-Fq sin - Fg sin 6g 


®C " “D 


Tilt Maneuver (Parking Lot Tilt) 

The maneuver modeled here is a modification of the tilt maneuver pre- 
sented in the previous memo. Here, is used as the value of angle- 
of-attack at a time half-way between the time to begin and end the tilt 
maneuver, Tj and Tp. Therefore, this procedure can fit any part of a 
parabola to the three points, depending on the value of and its rela 
tion to cq and ap. This is illustrated in Figure 2-2. 
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where otj^ 2 » “m3 three examples of values for which 
cause a different shape curve to be fitted. The calculation for 
during this maneuver is as follows: 

'V. 


“0 

T, 


OP 


r tan 


'BZ 

l'^bx j 


, calculated only at t = Tj 


= Tt * 


“d " “o 


(ar^ - a ) + (oLt-^ - a ) 

' D ' D 


lr]op 


“m + A,Z + A,Z‘ 
m 1 2 
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-N’p + 


m 


m 


O 0 


<^CG - 


yy 


No/m 


^2P ‘CGO * ® ■ '"26 


Zo - "3B 


BiKo " Bp ^pp 

ClP K3 - C^p Ki 


Aerodynamic Moment Control 


-^o Nq Pz 
hp ^'p 


Remaining Shaping and Control Parameters 

The following parameters are calculated after a value for Wq has been 
determined by some specified alpha policy. 
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= tan 



ec 


“r. - 6 

D R 


Vy 


'IP 


= -1 N' 
'iP P 

^yy 


' ^^0 ~ ^IP ^yy ' *' ^CGO ~ " "^CGB^ 






'DCG 


' ff ["o ^ "'p ■ ^0 'o ■ ‘bJ 


2.17.3 Nomenclature 


Vari able 

^1 ) ^2 
B > B 2 

C^p» C2p, C2g 

c 

^mo* 

^ZO’ ^Za 


Definition 

Temp, variables 
Temp, variables 
Temp, variables 

Mean aerodynanic cord 

Aero, moment coefficients 

Aero, normal force coefficients (pitch) 
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Variable 


Definition 


ax 


^2 * Sb 

•-0* *-B 

m 

No 

NV 


q 

s 

T 

T. 


'OP 


'^BX* '^BZ 
'^RX’ '^RZ 
^CG’ ^CG’ ^CG 

^Ig’ ^Ig 
^0’ ^0 
2 r 


DCG 


Drag moment due to CG offset 
Total force acting on vehicle 
Aeoo, axial force 
Moments of inertia 
Temp, variables 

Moment arms from CG to orbiter and booster engines 
Mass 

Normal force at zero angle-of-attack 
Partial of normal 
Dynamic pressure 
Aerodynamic reference area 
Present time 

Time to end tilt maneuver (begin alpha policy) 

Mid point between T^ and Tj 

Time to begin tilt maneuver (end vertical rise) 
Velocity components in body coordinates 
Inertial components of velocity relative to air 
Location of center of gravity 
Moment arm from CG to orbiter engines 
Average orbiter engine location 
Location of aerodynamic reference 
Location of body fixed accelerometer 
Temp, variable 

Desired body sensed acceleration 
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Variable 


Definition 


“D 


^CG8» ^CGO 


6 


D 



Desired angle-of-attack 
Value of angle-of-attack at Tgp 
Angle-of-attack at 

Deflection for thrust through center of gravity 
at booster and orbiter engines 

Desired orbiter engine deflection angle 

Gonmanded -pi tch;- angle- 

Acceleration at zero angle-of-attack 


2.17.4 Input/Output 

The inputs necessary for TSHAPE are the aerodynamic constants, configuration 
constants, I, m, CG’s, T, T^, T-j-, Vg, Vg, and Oj^. The outputs of this 
program are op, e^, «[,, Zj,,,. Sjug. and 
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2.18 TVC (Thrust Vector Control) 

2.18.1 Program Description 

This model describes the motions of massless engines with limits on 
deflection, deflection rate, and acceleration. 

2.18.2 Math Model 
Initialization (time = 0) 

<Sp “ 

foLD, 'nEH, 

''oLDi 'new, 
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Process and Limit Commands 


if (Sp - PNULL.) > PLIMIT. 

COMMAND^ . 


dL = PLIMIT. + PNULL. 

^COMMAND^ 


if (6„ - YNULL.) > YLIMIT. 

^COMMAND. ^ 


'COMMAND. = YLIMIT + YNULL. 


'^NEW, '^COHHAND, 


OLD 


I A K.1 rv I 


'NEW^ 


'COMMAND. 


OLD 


Calculate and Limit Rates 

= (5p - ^ )/Atinie 

^NEW^ ^NEW^ OLD^ 


limit to 6p 


LIMIT^ 


6v = (6. - Sy )/^time 

'newj new^ old,. 




limit to ‘’y 


LIMIT. 
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Cal culate and Limit Accelerations 


fip = (^b - \ 

^MEW. ^OLD. 


me 


limit to 6r 


LIMIT^ 


fiv - (K. • )/Atime 

''i ^NEW, ^OLD, 


limit to 6. 


'limit. 


Define new rate using limited acceleration 


6p = 6p + (L ) * (Atime) 
*^NEW^ '^OLD . ^NEW . . 


6y = 6y + (6y ) * (Atimo) 

^NEW^ ^OLD^. ^NEW. 


Define new position using limited rate 


6p =6 + (An ) * (Atime) + BIASP. 

^NEW. PqLD. *^NEW^ ’ 


6y = 6y + (6y ) * (Atlmc) + BIASY. 

’’new- ^OLD^. ^NEW^ 


90 


Reset Old Values 


6„ = <Sp - BIASP. 

^OLD. ’ 


= 5y * BIASY. 

Vd) old,. ’ 


^OLD. NEW^ 


«Y = <5y 

^OLD. ''new. 


Calculate Duty Cycle 


DCYCLE. = DCYCLE. + (Atime) * (l^p 1 + |5v D 
’ ^ *^NEW. ’'new. 

2.18.3 Nomenclature 
5 = present pitch deflection 


= present yaw deflection 


= past pitch deflection 


= past yaw deflection 


= present pitch deflection rate 


= present yaw deflection rate 



= past pitch deflection rate 



BIASP 

BIASY 

PLIMIT 

YLIMIT 

PNULL 

YNULL 

^P 

^COMMAND 

6y 

''command 

9 

aiMIT 


= past yaw deflection rate 
= pitch deflection acceleration 

= yaw deflection acceleration 
= engine pitch bias due to misalignments 

= engine yaw bias due to misalignments 

= pitch deflection limit 
= yaw deflection limit 
= pitch null position 
= yaw null position 

= pitch deflection command from flight computer 

= yaw deflection command from flight computer 
= pitch rate limit 


6w = yaw rate limit 

\lMIT 


^P 

\lMIT 


= pitch acceleration limit 
= yaw acceleration limit 
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i 


= engine index; 1, 2, 3 - orbiter engines 

4, 5 = solid engines 


6, 7 = OMS engines 

2.18.4 Input/Output 

Inputs: Pitch and yaw deflection commands from the flight control 

system (for each engine); deflection, rate, and acceleration 
limits; pitch and yaw null positions; and canned yaw bias. 

Outputs: Engine deflections, deflection rates, deflection accelerations, 

and duty cycle requirements. 



3.0 COORDINATE SYSTEMS 

Inertial Pblar-Equatorial - A right-handed orthogonal system with its origin 
at the center of the earth - X axis in the equatorial plane and positive 
through a reference meridian at the time of liftoff; the reference meri- 
dian is defined by the time of liftoff and the coordinate system used 
for gravity calculations. The Z axis is positive through the North 
Pole. 

Inertial Plumbline - An orthogonal system with its origin at the 
center of the earth, X axis parallel to the launch site gravity vector 
and positive in the direction opposite to gravitational acceleration. 

The Z axis lies in the launch plane and points downrange and the Y 
axis completes a right-handed triad. 

Local Vertical - An orthogonal system with its origin at the center 
of the earth, the X axis points from the earth center to the vehicle, 
the Z axis is in the plane containing the earth's rotation axis and 
the Xj^y axis. The Z axis is perpendicular to the X axis and points 
towards the North Pole. The Y axis completes a right-handed triad. 

Body ^ An orthogonal system with its origin at the engine gimbal 
pivot plane - X axis positive towards the nose of the vehicle along 
the main propellant tank centerline, Z axis positive "down", and 
the Y axis completes the right-handed system and is positive in the 
direction of the right wing. 
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Transformation matrix from polar-equatorial to plumbline coordinates 



11 

"12 

"13 

21 

®22 

"23 

31 

"32 

"33 


hi 

aj2 

®13 

®21 

®22 

®23 

"31 

®32 

"33 


= cos Xj_ cos 

/ * 

+ 

“e*L> 

* 

/ * 

“e*L> 

- cos-X|^ sin 



= sin 

~ sin sin X|^ cos (ugtj^+(f)|^) - cos sin (wgt|^+<{)|^ ) 

= sin A^^ sin A^^sin \ (t^gt[_+4>|_) 

* 

= -sin A|^ cos X|^ 

= -cos A|_ sin cos - sin A^^ sin (tOgt|_+<|>|_) 

“ic ^ ^ 

= -cos A|^ sin Aj_ sin (wgt|_+4i|_) + sin A|^ cos (ugt|^+<})|_) 
* 

= cos A|^ cos X|^ 


Where : 

X* = geodetic latitude of launch site 

(|)j^ = longitude of launch site 
Wg = angular rate of earth 
Tj^ = time of launch (from epoch) 

Aj^ s launch azimuth 
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Transformation matrix from body to inertial plumbline coordinates 



bjj = cos e cos Ip 

= sin e sin (p - cos e sin ip cos tj> 

b^3 = sin 8 cos <p + cos 6 sin ip sin <p 

bgj = sin Ip 

bg2 “ cos Ip cos <t> 

b23 = -cos Ip sin <p 

^31 ~ ® 


b 

b 


22 = cos 6 sin <p + 
22 = cos 0 cos <p - 


sin 9 sin ip cos <t> 
sin e sin ip sin (p 


where the Euler angles 8 , ip and are calculated in E 0 M. 
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Transformation matrix from local vertical to polar-equatorial coordinates 

[']■ 

djj = cos Xy cos <t’ 

di 2 = -sin 4' 

= -sin Xy cos 
d 2 j = cos Ay sin 
d 22 == cos <j>’ 

^23 “ 

d^^ sin Ay 

d32= 0 

^33 “ 

Where: 

sin Ay = Zp/R 
cos Ay =^1 - sin^ Ay 

sin 4>' = Yp/R cos Ay 

cos 4>* = Xp/R cos Ay 


11 

“12 

'^IS 

'21 

^22 

^*23 

bl 

'*32 

^*33 
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4.0 FLEXIBLE BODY MATH MODELS 

Mathematical models were developed for use in the SSFS. These models 
were programmed for a parallel burn solid rocket motor configuration. 
The flexible body version of SSFS is currently being checked out. 


4.1 Flexible Body Program Description 

This program contains the bending and slosh models for the launch 
configuration during first stage boost. It uses a generalized modal 
approach to bending which represents the elastic response by stan- 
dard normal modal equations with viscous damping. Included are models 
for aerodynamic forces and moment and thrust forces and moments to 
account for bending effects as well as the tail wags dog contribu- 
tion to bending. The rigid body and elastic response equations 
provided here are uncoupled and are considered separately since the 
magnitude of the coupling is insignificant. The number of equations 
is very sensitive to the vehicle configuration and to the complete- 
ness of the bending analysis. Therefore, when data becomes available 
it is likely that only a small percentage of the general set provided 
here will actually be required for SSV analysis. 

The model sums all the forces acting on each of the equivalent mass 
points and for a given mode numerically integrates the sum with a 
second order linear differential equation in modal displacement. 

The number of mass points at which aero forces and modal displacements 

are calculated will be less than 50. The number of modes at these 

points will be less than 10 each. The number of slosh masses will 

be less than 5 and the number of modes per slosh mass will be less 

than 5, The EOM, guidance atmosphere and control subroutine must be 

present to provide inputs for this model. 
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4.2 Vibration Equations 


N1 

Z (F .A.. + F .di.. + F - di.**) 
j = 1 ^ axj ^xij ayj azj 


+ Ml 

j ^^txj **’xij *^tyj '^yij ^tzj 

U 


+ K1 

Z (F .dh.. + F .d>**^F .d)..) 

J ' sxj ^XTJ syj ^yij szj ^zij^ 

J -A 


+ N2 

I (M .d>'..+M .<(>'• +M .d)'..) 


j = 1 


axj ^xij ayj "^yij azj ^zij 


+ M2 

r (M. . d>'-- + M. . d> ' • • + M. . it'--) 

j = 1 ^ 


+ K2 , • • 

Z (M •d>>«'^M • d«. • • ^ M . (t ^ • ) 

'' 'e v*i C\/T ^Vl T C71 ^7TT' 


j = 1 


sxj "^^xij syj ^yrj szj jzij' 


= m^. (q^. + 2 w, + u)? q.) 


q. = /q. dt + q^. 

q. = / q^. dt + q^. 
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Where: 


N1 

N2 

Ml 

M2 

K1 

K2 


0)* 




*^azj 

^txj 

^tyj 

^tzj 

^sxj 

'^syj 


■szj 

^axj 

^ayj 

"azj 


number of aero stations for aero forces 
number of aero stations for aero moments 
number of engines producing forces 
number of engines producing moments 
number of slosh stations for slosh forces 
number of slosh stations for slosh moments 
modal displacement due to bending mode i 
damping coefficient for mode i 
frequency of mode 1 
normalized mass for mode i 
aero forces in X direction at station j 
aero forces in Y direction at station j 
aero forces in Z direction at station J 
thrust forces in X direction for engine j 
thrust forces in Y direction for engine j 
thrust forces in Z direction for engine j 
slosh forces in X direction at station j 
slosh forces in Y direction at station j 
slosh forces in Z direction at station j 

aero moments about X -axis at station- j 

aero moments about y axis at station j 
aero moments about z axis at station j 
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= thrust moments about X axis due to engine j 

= thrust moments about Y axis due to engine j 

= thrust moments about Z axis due to engine j 

= slosh moments bbout X axis at station j 

= slosh moments about Y axis at station j 

= slosh moments about Z axis at station j 

= mode shape translation in X direction for mode i at location j 

= mode shape translation in Y direction for mode i at location j 

= mode shape translation in Z direction for mode i at location j 

= mode slope about X axis for mode i at location j 

= mode slope about Y axis for mode i at location j 

= mode slope about Zaxls for mode i at location j 

4.3 Aerodynamic Forces 


V 


0 

- wxpj. 

. - , - 


V 

" [a] [6] 

-''w ^w, 

J J 

-V,, cos A„, 

wypj 

V 

wzpj 


"j 


V' . . 
axbj 


» 

rv 1 
•^x 

Xp 


'^wxpj 

< 

V ' 

''aybj 



Ap 

- 

V 

wypj 


^azbj 

» « 



X 



. 

V 

wzpj 



M 


txj 

tyj 

tzj 


SXJ 

"szo 

♦xij 




^zij 

^iij 
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Vaxbj = Q (Zj ■ ^cg) - R 


V, u” = R (X,. - x^„) - P (z. - i) 
aybj ' j cg^ ' j eg 


»azb3 = P (''j ■ '^cg) - « - Xcg) 


.. M3 

''akbj ° j ♦xlj 'll 


V . '! = Z (t . . q . 
aybj j ^y^o 


'azbj = ,.^1 


'^axbj ^ '^aybj j ^zij ’ ^azbj Vi*j 


''axbj ^ qi 


^azbj “ '^axbj ^ '^yij “ '^aybj .^_ ^ '^xij 
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'^axbj “ '^axbj '^axbj '^axbj '^axbj 


I II i" IV 

V =V ,.+V L.+V .. + v .. 
aybj aybj aybj aybj aybj 


I II III IV 

'^azbj " ^azbj * '^azbj ^ '^azbj ^ '^azbj 


The previous 6 equations must be solved simultaneously for 


and The solution is as follows: 


let 


=‘11 “ 


‘12 


M3 

l 

i = 1 


‘•^zij ^i 


'13 


M3 

** * Q 4 

^ ^ 1 Vj 


M3 

®21 “ ** 


1 = 1 


ZIJ ^1 


®22 ■ 


M3 


^23 ~ *1 


1 = 1 


XU ^1 
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f 


M3 
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= Arctan ^ -y 


II 


- % P V,/ 

^axj 

= 5 <^XOj 

F . 
ayj 

■ ‘’j ^ 'yeo ®o 


F . = q. S (C . + C, . a.) 

azj ^0 '■ zoj zaj j' 


Where: 


S 

= magnitude of relative wind at Station j 

A 

= azimuth angle of relative wind at Station j 

'^wxpj ’ 


\yp j ’ 

\i 

= relative wind velocity at mass point j in platform 
coordinates 

W2pj 


[a] 

= transformation described in "coordinate systems" 

i 1 

= transformation described in "coordinate systems" 

[fi] 

= transformation described in "coordinate systems" 
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axbj* 
^aybj , 


'azbj 


= components of air 'Velocity in body coordinates 
at mass point j 


P, Q. R 




H 


V 

V 

V 


axbj’ 

aybj’ 

a 

azbj 


= angular velocity about X, Y, and Z axes respectively 
= location of mass point j in body coordinates 

= location of center of gravity in body coordinates 


= component of velocity at mass point j due to 
rotation of mass point about c.g. 


V “ 
^axbj * 


ir^ 


aybj’ 

V 

''azbj 


I 

'axbj’ 

I 

'aybj’ 

I 

azbj 


= components of velocity of vibrating mass with 
respect to rigid body 


= components of velocity due to perpendicular forces 
being rotated with respect to rigid body due to bending 


M3 = number of bending modes 

Vd »Vp ,Vp = velocity of vehicle relative to the earth in platform 
^Xp *^Yp ^Zp coordinate 

''axbj’ 

V . = components of velocity of mass point j with respect 

to air 

''azbj 
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’’j 


''aj 


P 

s 


= angle of attack of mass point j 
= mach number of mass point j 
= velocity of mass point j with respect to air 
= speed of sound 

= dynamic pressure at mass point j 

= mass density of air 
= aero reference area 




^Xaj ’ 




aero coefficients for mass point j 


sideslip angle for mass point j 
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cx> >< 


4.4 Engine Forces 


M3 




H3 


6 ‘ “ 0 * ♦ ♦ ^ 4 ^w*i •? 

pj pj i = 1 yij 1 


^bxa = eTj 


T. . = ET. sin 6„. 
byj J yj 


Tbzj ‘ -^3 *p3 V3 

Limit 0p^.j. and 6^^^ to position limit 


P3 =[“a=<Vcj - S3> * '^23 '^'p 3 ' «P3> * "l3 ^P3] 

Ki '[“a‘’<V3 ' *y3> * ^ ''13 ®'y3] /<'*"l3> 

Limit 0 to rate limit 

If e + edt exceed position limit, limit 0 to (position limit -8)/dt 

f 

®p3 ' ' %} ^ ®p3o 


e . = / 0 . dt + 0,,,.- 

yj yj yjo 


108 



M3 , 


ol 


0' . + Q + Z 'i’w! ,• Q-: 9^^■; 

pj ^ yy ep ep pj ep pj 


ep pj " ^^ecgj " 


0 


M3 


0* . + R + Z d»’ . . q. + 2 c . 0 ) 0’ . + 
yj 1=2 

„2 ; . + ^ (X . - X .) 

y^ ^2ze 


lU 0' . 

ey yj 


e' . = / 0 ' . dt + 0 ' . 
PJ ^ PJ PJO 


0' . = /e . dt + e' . 
PJ PJ PJO 


yj « * 


0' . = /0 . dt + e . 
yj yj yjo 


* ET. \ll + TAN2 e' + TAN^I 


txj J 
ftyj = 'tj ™ »'y 
ftzj’" -^tj ™N8'p 


e y 



M3 

^exj " “"’ej 

^ 1=1 

*^eyj " "'”ej 1 

[ \ " <*e3 


•■ M3 


XCG> - Ihj - ^CG> >" ^ 


M3 - -) 

i=j *yiJ ''1 J 
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. . , M3 .. 

^zj = ^ " ^cg^ ^ *zij '’i 


^txj ” ^txj '*’ ^exj 


^tyj " *"tyj ‘*‘ ”’ej " ^ecg^ ^ y* *^eyj 


*"tzj ■ ^tzj '’’ ^ej ^^ecg “ ^ p "*■ ^ezj 


Where: 



yaw engine gimbal angle with respect to the mounting 
surface for engine j 



= yaw engine gimbal angle with respect to rigid body 
coordinates for engine j 



pitch engine gimbal angle with respect to the mounting 
surface for engine j 



pitch engine gimbal angle with respect to rigid body 
coordinates for engine j 


location of engine j pivot point 
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location of engine j center of mass 



forces acting on vehicle 





= pitch engine actuator angle 
= yaw engine' actuator angle 

= pitch engine actuator rate 


I = moment of inertia of engine bell about Y axis at engine 

gimbal point 

I = moment of inertia of engine bell about Z axis at engine 

gimbal point 


6 . 
yj 


'pcj 

B . 

ycj 

^ep 


’ey 






= yaw engine actuator rate 

= pitch gimbal angle conmand 

= yaw gimbal angle command 
= damping factor for pitch engine dynamics 
= damping factor for yaw engine dynamics 
= frequency for engine dynamics 
= frequency for actuator dynamics 

= present engine rate for pitch 

= present engine rate for yaw 

= pitch engine angular acceleration 

= yaw engine angular acceleration 


= mass of engine j 



X 

z 



linear acceleration of vehicle in body coordinates 



= vehicle center of gravity 


The characteristicsfrequencies associated with the engine dynamics are 
much higher than the vehicle characteristic frequencies. It is 
recommended that the engine dynamics be integrated separately with an 
integration cycle of 50 milliseconds. 


4,5 Slosh Forces 

M3 .. • - 

“sj *xj ■*“ sj ^xj * J *xij 'll ■ \ 


* I* 


X. + 2c.co.X.+u 
yj ^sj “sj yj sj ^yj 


I '^yij ^^i " \ ^ ^^sj " ^cg^ 

i = 1 


+ P (Z^. - 

SJ eg' 
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2 « rt /V c- ^ 

‘zj 2 hi “sj ^zj * “ sj ^zj ' J +Z1J ll ■ ''z ■ ^'^Sj - ^cg> 


Q (X 


S 3 


- 

eg' 


^sxj “ '”’sj ^xj 


^xyj ‘ """sj ^yj 


^szj """sj ^2j 


^xj 


^'xj 


Vj 




ZJ 


^^ 2 j 


^Xj 


■^^Xj 


^yi 




^zj 


= /X 


ZJ 


Where: 

^xj» ^yj’ ^2j 
• • • 

^xj» Vj» ^zj 
*• »« #♦ 

^xj* Vj’ ^zj 


= displacements of slosh mass j 
= velocities of slosh mass j 
= acceleration of slosh mass j 
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= damping factor for slosh mode j 
= characteristic frequency of slosh mode j 


position of slosh mass j 


m . = mass of sloshing fluid at, mode j 


4.6 Aerodynamic Moments 


^axj ' ^ayj ^^cg " ^arj^ " .^azj ^^cg " ^arj^ 


^ayj ■ ^azj ^^cg " ^arj^ " *^axj ^^cg “ ^arj^ 


''azj = %xj - ''arj) ‘ %} (><cg ' ><arj) 


«axj - »> «:i6j ) 


^ciyj ==Q-S C (C . + C .a.> 
^ nioj maj J > 
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M . = q. S b C • 
azj n 3 J 


^axj ' ^axj + 


%j ' ”ayi + 


M . + M I + M " 
azj azj azj 


Where; 


^l 6 a* ^nsa’ 


= aero coefficients for station j 


6 - = aileron deflection 

a 


4.7 Engine Moments 


^txj "^ybj ^^ej " ^zbj 


“^zbj ^^ej ' ^cg^ ^xbj ^^ej " ^cg^ 


= -^xbj <^ej - ''eg) * Vj <*ed ' ><09) 
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^txj 


M3 M3 

'"^ybj f 1 ^zij ‘’i ^ ^ j Vij ^i 


M '! 

■^tyj = -T 


M3 ^ + T M3 . „ 

zbj I "^xid ^xbj z '^zij 
i = 1 i = 1 


M. ■! 
tzj 


_T M3 
" xb j z 


d, . . q. + T 
^yij ..u.- 


1 = 1 


M3 


ybj z 

i = 1 


'’i 


^^exj ~ ^ *^ezj'”’ej^^ej“^ecg^*©P0 ^ ^^ej'^cg^“^^eyj''^ej^^ej'^ecgo^ *®yj^^^ej”^cg^ 
^ey j ' ^ex j ^ ^e j ‘^ecg j ^ ^ ^ez j '"'e j ^ ^e j “^ecg j ^ p j ^ ^ ^e j '^cg ^ 


*^ez j ^ ^ey j ^e j ^ ^e j "^ecg j ^ y j ^ j "^cg ^ '^ex j ^ ^e J '^cg ^ 


*^txj ~ *^txj ^txj ^ ^exj 


M* . = M* '. + M* " + M . 

tyj tyj tyj eyj 


^tzj = "tzj + $ M^^ . 


Where: 


I . = moment of inertia of engine j 


116 



4.8 Slosh Moments 


'^sxj = \ *20 - ''z ™sj *yj 


M 

syj 


'"sj ^xj 




M I 
szj 


^ ^yj 


\ ^a 'xi 


*^sxj ~ *^syj “ ^sj^ “ ^szj 


^syj “ *'szj ^^cg ” " *^sxj ^^cg " ^sj^ 


«sz3 = ^xa (^cg - Y 3 J) “ ^yj ^^cg ‘ ^sj) 


M = M ' + M '.' 
"sxj ^SXJ 


^syj ■ '^syj * ^syj 


M . = M '. + M '1 

SZJ SZJ SZJ 


117 



